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Abstract 



A number of positive and null results on the time variation of fundamental 
constants have been reported. It is difficult to judge whether or not these 
claims are mutually consistent, since the observable quantities depend on 
several parameters, namely the coupling strengths and masses of parti- 
cles. The evolution of these coupling-parameters over cosmological history 
is also a priori unknown. A direct comparison requires a relation between 
the couplings. We explore several distinct scenarios based on unification of 
gauge couplings, providing a representative (though not exhaustive) sam- 
' pie of such relations. For each scenario we obtain a characteristic time de- 

\ pendence and discuss whether a monotonic time evolution is allowed. For 

t^- ■ all scenarios, some contradictions between different observations appear. 

We show how a clear observational determination of non-zero variations 
OO ' would test the dominant mechanism of varying couplings within unified 

\ theories 

OO 

o 



1 Introduction 



Any observation of a time variation of "fundamental constants" would be a far- 
reaching discovery. There are various claims for a detection, and many more 
observations indicating a null variation: criteria for judging their mutual consis- 
tency would be useful. We investigate whether a simple scenario exists which 
can account for several observational claims simultaneously and is consistent 
with unification of Standard Model (SM) gauge couplings. Several observations 
motivate such a study. First, the claimed deviations from the present value of 
the fine structure constant a, or the proton-to-electron mass ratio fj,, observed 
in quasar absorption systems. Second, the discrepancy between the primordial 
7 Li abundance expected from standard nucleosynthesis (BBN) and seen in old 
halo stars, which may be explained by a variation of "constants" within a uni- 
fied framework [HE]- Third, the theoretical insight that scalar fields cannot 
be exactly constant over the entire cosmological evolution, and that a possible 
"late" time evolution can play an important role in the dynamics of the ex- 
panding Universe. A time variation of couplings arising from the evolution of 
a "cosmon" field in so-called "quintessence scenarios" O H] would link these 
variations to observables in cosmology [51 E] . 
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Due to the many unknowns of the underlying particle physics models and 
the partly contradictory present observational situation, a systematic treatment 
is not easy, and may even seem premature. Nevertheless we consider a first 
attempt to be useful, in order to discuss strategies that can be used to compare 
variations of different observables. The power of the proposed method will only 
become clear if and when future observations present a less ambiguous picture. 

The basic approach in the present paper relates the fractional variations of 
different fundamental couplings G^, such as the fine structure constant a, the 
proton-electron mass ratio \i or the ratio of the nucleon mass to the Planck mass, 
by an assumption of proportionality, with fixed "unification coefficients" d^. 
The choice of the values of dk is in turn determined within different scenarios of 
varying parameters in unified theories (GUTs) where the gauge couplings of the 
Standard Model converge at a unification scale Mx- The assumption of time- 
independent coefficients dk covers a large class of possible models for varying 
couplings. This assumption is, however, not a necessity, and we will describe 
specific quintessence models where it may not be realized in a forthcoming 
paper [TJ. 

In Section [2] of this paper we review observational determinations of the 
variation or constancy of couplings, considering five types of methods: early- 
universe cosmology, astrophysical spectroscopy, nuclear physics in the Earth and 
the Solar System, gravitational physics, and atomic clock comparisons. In Sec- 
tion [3] we introduce the unification of couplings (Grand Unified Theory, GUT) 
and determine the implications of unification and supersymmetry (SUSY) for 
the Standard Model parameters and for the observables we consider. We further 
define six unified scenarios by considering different possibilities for the variation 
of the Fermi scale and the superpartner masses. Within each scenario we reduce 
the various observational results to constraints on the time evolution of a single 
fundamental coupling, and discuss the mutual consistency of observations. 

In Section H] six cosmological epochs are introduced, and the constraints on 
variation deduced in Section [3] are collected into a set of "evolution factors" 
for each unified scenario. These evolution factors are a measure of the overall 
size of coupling evolution between a given epoch and the present. We then 
determine for each scenario to what extent a monotonic variation over time can 
be consistent with the data. Section [5] draws some general conclusions. 

In a subsequent paper [7] we will investigate the scalar field dynamics that 
could give rise to a small but nonzero variation of couplings. The presence of 
a cosmologically varying degree of freedom gives rise to important additional 
effects. It affects gravity on large scales, altering the expansion of the Universe 
and potentially giving rise to the observed late-time acceleration. Also on local 
scales, a light field weakly coupled to matter produces long-range forces which 
are tightly constrained [8] by Solar System precision tests of gravity and the 
null results of experiments testing the Weak Equivalence Principle. Combining 
all these considerations leads to tighter constraints on models but also offers 
more possibilities to test them. 
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2 Data: variations and constraints 



Here we review and discuss the observational data that we will consider in our 
effort to obtain a unified picture of time variation of couplings. We summarize 
the results that are most relevant for our analysis in Table [TJ 

2.1 Early universe: BBN and CMB 

The earliest processes for which Standard Model physics can be tested are BBN 
(z ~ 10 10 ) and CMB (z ~ 10 3 ). Hence they constitute the most far-reaching 
tests of a possible variation of couplings. 

BBN The influence of varying constants on BBN has been studied extensively 
[21 El [lOl [III [121 [131 [III [IS [111 [171 [111 [19] (see also [20] and references therein) . 
The method developed in [TJ [12] accounts for possible simultaneous variations 
of different fundamental parameters, while previous studies restricted attention 
to the variation along one particular direction in parameter space. 

Our approach in pQ was first to calculate the leading dependence of BBN 
abundances on a set of "nuclear parameters", comprising elementary particle 
coupling strengths and masses and nuclear binding energies. In a second step, 
these nuclear parameters were related to fundamental parameters in particle 
theory, which allowed us to consider (at linear order) any combination of vari- 
ations at BBN. Thus, our results are independent of any assumptions about 
unification. 



In an extension of our previous treatment, we include in Appendix A the 
possible effect of varying constants at CMB on the input parameter rj of our 
BBN procedure. The important parameter is here the variation of miy/Mp at 
CMB relative to the variation at BBN, where is the nucleon mass and Mp 
is the "reduced" Planck mass. We consider two limiting cases. First, when 
A(mjv/Mp)| CMB <C A(m7v/Mp)| BBN : then our previous results hold. In the 
second case, with A(totv/Mp)| CMB ~ A(mAr/Mp)| BBN , the value of r\ may be 
significantly rescaled. 

Our result for the leading dependence of primordial abundances Y a on fun- 
damental particle physics parameters G/% = (G,a, {4>},ni e ,5 q ,rh) is summa- 
rized in Table [2j Here 5 q = rrid — m u denotes the light quark mass difference, 
rh = (rrid+Tn u )/2 the average light quark mass, and {(ft) is the expectation value 
of the Higgs scalar that determines the Fermi scale of the weak interactions. 
The variation of dimensionful quantities is defined relative to the QCD strong 
coupling scale A c . We found that only D, 4 He and 7 Li can be used to constrain 
parameters at BBN. Whilst Table [2] only gives linear dependencies, we can 
account for nonlinearities by running the full BBN code with the appropriate 
variations pQ. 

The current observational and theoretical values for the Y a are given in 
Table EJ where Y p is the helium mass fraction equal to four times the ratio 
of 4 He number density to hydrogen. The uncertainty in the rj determination, 
n = (6.20 ± 0.16) x 10" 10 (WMAP5 plus BAO and SN, [21]) yields a further 
correlated error for the abundances, which can be treated using the method of 
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Method 


redshift 


A In a 

[io- 6 ] 


[io- 5 ] 


A In Gm 2 N 

[io- 2 ] 


Alnx 

[io- 5 ] 


Alny 

[io- 5 ] 


AlnF 

[io- 5 ] 


A In F' 

[io- 4 ] 


AlnAi 87 
[IO" 2 ] 


Oklo a [45J 


0.14 


0.00 ± 0.06 
















21cm [H] 
Sun [57] 
Heavy/Hi, low-z 02] 
Meteorite 09] 
Ma epoch 2 [28] 
Ammonia [39j 

21cm IH] 
HI / OH [43] 


0.247 
0.43 
0.40 
0.44 
0.65 
0.68 
0.685 
0.765 


-2.9 ±3.1 


0.06 ±0.19 


0±0.72 


1.0 ± 1.7 


-0.20 ±0.44 
-0.16 ±0.54 


0.4 ± 1.1 




3.3 ±3.2 


Absorption [31] 
Ma epoch 3 [28] 
Absorption [32J 
Heavy/Hi, high-z 02] 


1.15 
1.47 
1.84 
2.03 


-0.1 ± 1.8 
-5.8 ±1.3 
5.7 ± 2.7 






0.6 ± 1.9 










#2 [37] 
Ma epoch 4 [28] 

#2 E7J 
Neutron stars [60] 

en / co bu ' 

CII / CO [44] 


2.59 
2.84 
3.02 
3.3 
4.69 
6.42 


-8.7 ±3.7 


2.78 ±0.88 
2.06 ± 0.79 


-0.7 ±2.4 








1.4 ± 1.5 
0.1 ± 1.0 




CMB [25, [26] 


10 3 


n +ixio 4 

-3xl0 4 




0+ 7 













Table 1: Observational la bounds on variations. Observables are defined as \i = m p /m e , x = a 2 g p fi , y = a 2 g Pl F = g p [a 2 fj] , 
F' = a 2 1 [i. The given redshift may denote a single measurement, or an averaged value over a certain range: see main text. The two 
CMB bounds are independent of each other. Our BBN bounds cannot be displayed in this form. 





D 


4 He 


7 Li 


G 


0.94 


0.36 


-0.72 


a 


3.6 


1.9 


-11 


{<t>) 


1.6 


2.9 


1.7 


m e 


0.46 


0.40 


-0.17 


s g 


-2.9 


-5.1 


-2.9 


m 


17 


-2.7 


-61 


7] 


-1.6 


0.04 


2.1 



Table 2: Sensitivity of abundances Yi to variations of fundamental parameters 
Gk and the baryon-to-photon ratio r\. 



Abundance 


Observational 


Theoretical 


D/H 

Y P 
7 Li/ H 


(2.8 ±0.4) x 10~ 5 

0.249 ± 0.009 
(1.5 ±0.5) x HT 10 


(2.61 ±0.04) x 10~ b 

0.2478 ± 0.0002 
(4.5 ±0.4) x 10" 10 



Table 3: Current observational and theoretical primordial abundances 



22]. For any given set of fundamental variations we can define 

^(Yi-Yf^w^-Y? 



X' 



'■J 



with the inverse weight matrix 



w. 



a 



2.,, 



+ 5. 



/_ 2 
l] \ a obs,i 



-1 



where 



We take, as in pQ, 



a 



2,V 



YY d\nYid]nYj 
1 3 dliir] d\ni] 



dln(D/H,y p , 7 Li/H) 



(-1.6, 0.04, 2.1). 



(1) 



(2) 



(3) 



(4) 



din 7/ 

The 1(2)<7 error contour is given by x 2 /^ — 1(4) where v is the number of 
degrees of freedom. As the final abundances depend on variations of all funda- 
mental constants, we have to evaluate the variations allowed by BBN for every 
model separately. 

It has been pointed out that the important 8 Be resonance very near the 
ground states of 7 Be±n and 7 Li±j> makes the exchange reaction converting 7 Be 
into 7 Li potentially sensitive to variations in nuclear forces. We give an estimate 



of this sensitivity in Appendix B and show that it is unlikely to be significant 
for the range of variations that we consider. 

In the light of complex astrophysics which may affect the extraction of 
the primordial Li fraction [23] , we also consider bounding the variations using 
deuterium and 4 He alone. This yields a value consistent with zero for variations 
at BBN, since these abundances are consistent with standard BBN. 
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CMB In principle, a and G are bounded by CMB, but there are significant 
degeneracies with other cosmological parameters [24, 25J: see also the discussion 
in Appendix A. Current bounds are 

0.95 < < 1.02 (2a). (5) 

The CMB anisotropies may also be used to constrain the variation of Newton's 
constant G. The resulting bound depends on the form of the variation of G 
from the time of CMB decoupling to now. Using a step function one finds 
[261 [27] 

G 

0.95 < — < 1.05 (2a), (6) 
Go 

where the instantaneous change in G may happen at any time between now 
and CMB decoupling. Using instead a linear function of the scale factor a, the 
bound is 

G 

0.89 < — < 1.13 (2a). (7) 
Go 

Note that here, as in most studies of time-dependent G, units are implicitly 
defined such that the elementary particle masses (and thus the mass of gravitat- 
ing bodies, if gravitational self-energy is neglected) are constant. The relevant 
bound on dimensionless parameters concerns Gm 2 N = (m^ / M-p) 2 (&ir)~ l . 

2.2 Quasar absorption spectra 

The observation of absorption spectra of distant interstellar clouds allows to 
probe atomic physics over large time scales. Comparing observed spectra with 
the spectra observed in the laboratory gives bounds on the possible variation 
of couplings. Different kinds of spectra (atomic, molecular, . . . ) are sensitive 
to different parameters, principally a and \i = m p /m e . 

Atomic spectra are primarily sensitive to a. Several groups using various 
methods of modelling and numerical analysis have published results; we quote 
here only the latest bounds. Murphy and collaborators [28j studied the spectra 
of 143 quasar absorption systems over the redshift range 0.2 < z a f, s < 4.2. Their 
most robust estimate is a weighted mean 

^ = (-0.57 ±0.11) x 10~ 5 . (8) 

In discussing unified models in Section 01 we will define various "epochs" 
for the purpose of collating data and comparing them with models over certain 
ranges of redshift. The 143 data points are then assigned to different epochs: 
we choose to put boundaries at z = 0.81 and z = 2.4, thus we obtain three 
sub-samples 

z < 0.81, N svs = 18, (z) = 0.65, — = (-0.29 ± 0.31) x 10~ 5 

a 

0.81 < z < 2.4 N sys = 85, (z) = 1.47, = (-0.58 ± 0.13) x 10~ 5 

z > 2.4, N sys = 40, (z) = 2.84, = (-0.87 ± 0.37) x 10~ 5 . (9) 
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Here we have used the "fiducial sample" of [29], the weighted average has been 
taken, and we have included [30] the 15 additional samples used in [28] . For 
convenience we will refer to these results as "Ma". 

Further results have been obtained by Levshakov et al. [32], and reported 
in [31]: 

— = (-0.01 ±0.18) x 10" 5 , z abs = 1.15 
a 

— = (0.57 ± 0.27) x 10 -5 , z abs = 1.84. (10) 
a 

Previously [33\ [34"] more stringent null results were claimed, but doubts have 
been cast [35] (see also [36]) on the validity of the statistical analysis. We 
note that the value for z = 1.84 has an opposite sign of variation to the Ma 
result, though the variation does not have high statistical significance. The 
observational situation is clearly unsatisfactory. 

Vibro-rotational transitions of molecular hydrogen H% are sensitive to \i = 
m p /m e . From H2 lines of two quasar absorption systems (at z = 2.59 and 
z = 3.02) a variation is found [37] of 

^ = (2.4 ±0.6) x HT 5 , (11) 
A* 

taking a weighted average. We will refer to this result as "R/i" after Reinhold 
et al. The individual systems yield [37] 

— = (2.78 ± 0.88) x 10" 5 , z abs = 2.59 
A* 

— = (2.06 ± 0.79) x 10" 5 , z abs = 3.02. (12) 
A* 

Recently the z = 3.02 system has been reanalyzed |38j . with the result that the 
claimed significance of Eq. (|12p was not reproduced, and the absolute magnitude 
of the variation is bounded by |A/i//z| < 4.9 x 10~ 5 at 2a, or 

\A(i//ji\ < 2.5 x 10~ 5 , z abs = 3.02 (la). (13) 

The inversion spectrum of ammonia has been used to bound /i precisely at 
lower redshift [39]. Recently the single known NH3 absorber system at cosmo- 
logical redshift has been analysed [40], yielding 

— = (0.74 ±0.89) x 10~ 6 , z = 0.68. (14) 
A* 

This is a considerably stricter bound but applies at a different epoch. Extrapo- 
lation to today with linear time dependence gives fi/fi = (1.2± 1.4) x 10 -16 y^ 1 . 

The 21cm HI line and molecular rotation spectra are sensitive to y = a 2 g p , 
where g p is the proton g- factor. Bounds on this quantity from [41] are 

— = (-0.20 ±0.44) x 10~ 5 , z = 0.247 

y 

^ = (-0.16 ±0.54) x 10" 5 , z = 0.685. (15) 

y 
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Further, the comparison of UV heavy element transitions with HI line probes 
for variations of x = a 2 g p /j,~ 1 [42J: the weighted mean of nine analysed systems 
yields 

Ax 

— = (0.63 ± 0.99) x 10~ 5 , 0.23 < z abs < 2.35. (16) 

However, we note that i) the systems lie in two widely-separated low-redshift 
(0.23 < z < 0.53) and high-redshift (1.7 < z < 2.35) ranges; and ii) these 
two sub-samples have completely different scatter, x 2 /^ about the mean for 
the low- and high-redshift systems being 0.33, and 2.1, respectively. Hence we 
consider two samples, with average redshift z = 0.40 (5 systems) and z = 2.03 (4 
systems). With expanded error bars in the high-redshift sample (after "method 
3" of 02]) we find 

— = (1.02 ± 1.68) x 10~ 5 , (z) = 0.40 

x 

— = (0.58 ± 1.94) x 10" 5 , (z) = 2.03. (17) 

x 

The comparison of HI and OH lines is sensitive to changes in F = g p [a 2 fi] 157 
[43J and yields 

A F 

— = (0.44 ± 0.36 stat ± 1.0 sys ) x 10~ 5 , z = 0.765. (18) 

A similar method comparing CII and CO lines has very recently been proposed 
at high redshift [H] yielding the best bound at redshifts > 4.5. The following 
bounds on F 1 = a 2 //i are obtained for two systems: 

AF' 

-yj- = (0.1 ± 1.0) x 10" 4 , z = 6.42 
Ai* 1 ' 

-yj- = (1.4 ± 1.5) x 10" 4 , z = 4.69. (19) 

2.3 Terrestrial and Solar System nuclear constraints 

Oklo natural reactor From modelling nuclear reaction processes which hap- 
pened in the Oklo mine about two billion years ago (At ~ 1.8 x 10 9 y, z ~ 0.14 
with WMAP5 best fit cosmology) one can in principle bound the variation of 
a over this period. The determination of A In a at the time of the reactions 
results from considering the possible shift, due to variation of electromagnetic 
self-energy, in the position of a very low-lying neutron capture resonance of 
149 Sm. The analysis of [45] gives the bound (taken as la) 

- 5.6 x 10~ 8 < Aa/a < 6.6 x 10~ 8 . (20) 

For a linear time dependence this results in the bound 

\a/a\ < 3 x 10" 17 y -1 - 

A more recent analysis using different reactor models and consequently different 
neutron spectra [46] found 

-2.4 x 10~ 8 < Aa/a < 1.1 x 10~ 8 
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with an additional nonzero solution (due to the other branch of the resonance 
peak) at Aa/a ~ 8 x 10~ 8 . We will consider the more conservative bound. 

Note that these results concern varying a only. If other parameters affecting 
nuclear forces, in particular light quark masses, are allowed to vary, the inter- 
pretation of this bound becomes unclear [47\ [48J since it depends on a nuclear 
resonance of 150 Sm whose properties are very difficult to investigate from first 
principles. In the absence of a resolution to this problem we consider Oklo 
as applying only to the a variation in each model. In scenarios where several 
couplings vary simultaneously we do not consider strong cancellations. Never- 
theless, we allow for a certain degree of accidental cancellation and therefore 
multiply the error on the bound Eq. (|20D by a factor three. 



Meteorite dating Long-lived a- or /3-decay isotopes may be sensitive probes 
of cosmological variation [49|. Wl\ [50] . Their (generally) small Q- values result 
from accidental cancellations between different contributions to nuclear bind- 
ing energy, depending on fundamental couplings in different ways, thus the 
sensitivity of the decay rate may be enhanced by orders of magnitude. 

The best bound concerns the 187 Re /3-decay to osmium with Qp = 2.66 keV. 
The decay rate Ai87 is measured at present in the laboratory, and also deduced 
by isotopic analysis of meteorites formed about the same time as the solar 
system, tyiet — 4.6 x 10 9 y ago (z ~ 0.44). More precisely, the ratio Ai87/Au, 
averaged over the time between formation and the present, is measurable |49} 
[51], where Au is the rate of some other decay (for example uranium) used to 
calibrate meteorite ages. The experimental values of Ai87 imply (setting Xjj to 
a constant value) 

*mL / A \ 187 ^ dt = 0.016 ± 0.016. (21) 

J-t Mct Ai 87 

Since the redshift back to tyiet is relatively small, we obtain bounds on recent 
time variation by assuming a linear evolution up to the present, for which the 
LHS is — (tMet/2)Ai87/Ai87 and the fractional rate of change is bounded by 

Al87 , _ „ , „ ^ 12 



Al87 



(-7.2 ±6.9) x lO-^y- 1 . (22) 



Projected back to tMet this gives the bound A In Ai87 — 0.033 ±0.032 (z ~ 0.44). 
This is a conservative bound unless the time variation has recently accelerated 
(which we consider unlikely), or there are significant oscillatory variations over 
time. 

The decay rate varies as 



Ai 87 oc G\Q%ml oc (</>) 2 Q%y 2 e 

where y e is the electron Yukawa coupling. Variation in Qp is determined by 
the near-cancellation between the nuclear Coulomb self-energy and asymmetry 
energy of 187 Re and 187 Os, and the nucleon masses, via 

AQ/3 ~ 0.77Aa A - 26Aa c + A(5 N - m e ), 
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where a a — 23.7 MeV and ac — 0.71 MeV are the asymmetry and Coulomb 
terms of the semi-empirical mass formula and 5n = m n — m p is the nucleon 
mass splitting. Thus the fractional variation in A187 is 

Aln ^1-2.1 x l0 4 Aln^-2.1 x 10 4 A In a + 880 A In * N - me . (23) 

Since the possible dependence of "control" decay rates \\j/mN on nuclear or 
fundamental parameters is much weaker than that of A187 / m K? we use this 
result for the variation of A187 in units where Au is constant! 1 ] Then using 
relations derived in [52] and considering also the effect of varying m s on the 
nucleon mass, 

A In JjV ~ me ~ 2.6 A In ^ - 0.65 A In ^ - 0.97 A In a - 0.12A In ^ , (24) 
m N A c A c A c 

A In « -0.9 A In H (25) 
m N A c 

we find the decay rate dependence to be 

A In — ~ -2.2 x 10 4 Alna- 1.9 x 10 4 A In — + 2300 A In - 580 A In — ^. 

m N A c A c A c 

(26) 

2.4 Bounds on the variation of the gravitational constant 

Variations of Newton's constant have been studied in the solar system and in 
astrophysical effects. Whilst all references give bounds exclusively on a potential 
variation of G, one should note that besides G also nuclear parameters (neutron 
/ proton masses and parameters of nuclear forces) can vary, which would in 
general add degeneracies and make the results less stringent. It has generally 
been assumed that particle masses are constant, thus the resulting bounds 
actually constrain variation of Gm 2 N oc (m^v/Mp) 2 . 

In the solar system, changes of G induce changes in the orbits of planets. 
Range measurements to Mars from 1976 to 1982 can be used to obtain [53] 



G/G = (2 ±4) x 10- 12 y _1 - (27) 

Lunar laser ranging from 1970 to 2004 yields [54j 

G/G = (4 ±9) x 10- 13 y _1 . (28) 

The stability of the orbital period of the binary pulsar PSR 1913+16 [55] may 
be used to deduce 

G/G = (1.0 ± 2.3) x 10~ n y -1 . (29) 

Recent observational advances may improve such bounds considerably, with a 
limit of 

G/G = (0.5 ± 2.6) x 10~ 12 y _1 (30) 



1 Variation of A c alone would not cause a dominant effect on A187. Both the Coulomb and 
asymmetry terms scale with A c , thus the effect of varying A c is confined to the last term on 
the RHS of (|23|) . i.e. varying the ratio of the weak scale to the strong scale. 



10 



from PSR J04374715 quoted in the preprint [56]. All these results apply at the 
present epoch z = 0. 

A bound on the behaviour of G over the lifetime of the Sun (approximately 
4.5 x 10 9 y, z = 0.43) was found by Krauss et al. [57] by considering the effect of 
the resulting discrepancy in the helium/hydrogen fraction on p-mode oscillation 
spectra. The claimed constraint is 

\G/G\ < 1.6 x Hr 12 y -1 
|AlnG| < 7.2 x 1(T 3 z = 0.43, (31) 

where the assumed form of variation is a power law in time since the Big Bang, 
which may be approximated over the last few billion years as a linear depen- 
dence. For models with significantly nonlinear time dependence the bound may 
be reevaluated: since the bound arises from the accumulated effect of hydrogen 
burning since the birth of the Sun, it may be expressed as an integral of the 
variation over the Sun's lifetime analogous to Eq. (|2ip . 

The properties of compact objects such as white dwarfs and neutron stars 
(NS) have been used to bound possible variations of G: see for example [58jl59]. 
The strongest bound not relying on speculative physical effects arises from 
comparing the masses of young and old neutron stars in binary systems |BD] : if 
one member of the binary is a pulsar, precision timing can be used to determine 
the masses. The mass of neutron stars at formation is determined to first 
approximation by the Chandrasekhar mass 

Mc - * chi (32) 

where m n is the neutron mass. This may be reexpressed in terms of the baryon 
number of the star rig cx Mch/™>n oc (Gm 2 ) -3 / 2 , which is expected to be con- 
served up to small corrections from matter accreting onto it. Thus the relative 
masses of NS measured at the same epoch probes the fractional variation of 
Gm 2 between their epochs of formation. Given that the oldest neutron stars 
are up to 12 Gy old, z ~ 3.3, the variation of the average NS mass \i n is found 
to be (i n = -1.2 ± 4.0(8.5) x l^ l2 M & y- 1 at 60% (95%) confidence level. In 
units where particle masses are constant, we have 

G/G = -0.6 ± 2.0 (4.2) x 10~ 12 y -1 , (33) 

where the averaging is performed over the last 12 x 10 9 y, and the bound should 
be reinterpreted for variations which are not linear in time. The absolute vari- 
ation over this period is then bounded at ler as 

A In G = (-0.7 ± 2.4) x 10~ 2 , z = 3.3. (34) 
2.5 Atomic clocks 

Stringent bounds on the present time variation of the fine structure constant and 
the electron-proton mass ratio may be obtained by comparing different atomic 



11 



transitions over periods of years in the laboratory [61]. A recent evaluation [62] 
of atomic clock data gives 

dhxa/dt = (-0.31 ± 0.3) x 10~ 15 y _1 

dhin/dt = (1.5 ±1.7) x Hr 15 y -1 - (35) 

Fortier et al. [63j obtain stronger bounds, \a/a\ < 1.3 x 10~ 16 y _1 , if other 
relevant parameters are assumed not to vary. If other atomic physics parameters 
are allowed to vary, this bound becomes considerably weaker, depending on a 
possible relative variation of the Cs magnetic moment and the Bohr magneton. 
Direct comparison of optical frequencies may yield bounds at the level of 10~ 
per year; limits on variation of a from this method are reported with uncertainty 
2.3 x 10~ 17 /y |64j but designated as preliminary. If these bounds are used then 
our limits from atomic clocks via a variation should be tightened by about an 
order of magnitude. 

Extrapolating the results of [62J to the time of Oklo [z = 0.14, t = 1.8 x 
10 9 y) gives 

Alna = (-0.56 ±0.54) x 10~ 6 , 

Aln/i = (-0.27 ±0.31) x 10" 5 . (36) 

The bound on a at this epoch is weaker than that from Oklo, Eq. (|20p . The 
bound on [i cannot be directly compared, due to unquantified theoretical un- 
certainty in the Oklo bound. However, if we interpret Oklo as simply bounding 
Aa/a, we find that it provides the strongest bound on fi variation for all uni- 
fied scenarios we consider (see Section [3} except our scenario 3, where the high 
ratio Aln/i/Alna ~ —325 means that atomic clock bounds on \i are the most 
sensitive. 



3 Unification and relations between 
coupling variations 

In this paper we consider the hypothesis that, for all redshifts, all fractional 
variations in the "fundamental" parameters (see section 12. 1|) are propor- 
tional to one nontrivial variation with fixed constants of proportionality. If the 
variation of the unified gauge coupling A In ax is nonvanishing, we may write 

AlnGfe = d k A\na x (37) 

for some constants dk, assuming small variations. Different unification scenar- 
ios correspond to different sets of values for the "unification coefficients" d k . 
Considering the values of A In (7/; as coordinates in an iV^-dimensional space, 
this assumption restricts variations to a single line passing through zero. The 
variation then constitutes exactly one degree of freedom. We will go beyond 
this hypothesis in a subsequent paper [7j where we consider a specific model for 
which a fixed linear relation (|37p is not realized for all z. 
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3.1 GUT relations 



It is natural in this context to consider models with unification of gauge cou- 
plings (GUTs). These have the property that variations of the Standard Model 
gauge couplings and mass ratios can be determined in terms of a smaller set 
of parameters describing the unified theory and its symmetry breaking. Hence, 
if nonzero variations in different observables are measured at similar redshifts, 
models of unification may be tested without referring to any specific hypothesis 
for the overall cosmological history of the variation. We need only assume that 
for a given range of z the time variation is slow and approximately homogeneous 
in space, hence A In ax depends only on redshift z to a good approximation. 
The relevant unified parameters are the unification mass Mx (relative to the 
Planck mass) , the unified gauge coupling ax defined at the scale Mx , the Higgs 
v.e.v. ((f)) and, for supersymmetric theories, the soft supersymmetry breaking 
masses rh, which enter in the renormalization group (RG) equations for the 
running couplings. Then, for the variations at any given z we can write 

A\n^ = d M l, Alna x = d x l, Aln^- = d H l, Aln^- = d s l, (38) 

where l{z) is the "evolution factor" introduced for later convenience. If ax 
varies nontrivially we may normalise I via dx = 1. In supersymmetric theories 
we set ax = 1/24, in nonsupersymmetric theories we set ax = 1/40 and ds = 
[65]. 

We make the simplifying assumption that the masses of Standard Model 
fermions all vary as the Higgs v.e.v., i.e. Yukawa couplings are constant at the 
unification scale: 

Aln^ = Aln-^ = Aln-^ = Aln^ = Aln#. (39) 
M x M x M x M x M x 

Here we neglect the effects induced by a variation of ax on the RG running 
of fermion masses, and consider the quark masses defined at an appropriate 
RG scale for low-energy observables. We have explicitly calculated the effect 
of varying a%{Mx) on the running of quark masses: for low-energy observables 
such as m q (Q 2 ) / A c one should consider an RG scale Q 2 that is fixed relative to 
A c . Thus the variation of m q (Q 2 )/m q (M x ) is entirely due to the dependence 
on as(Mx), which is suppressed by a loop factor ajc/i compared to the non- 
perturbative dependence of A c /Mx on axil Such effects enter at the order of 
1% correction, which is already smaller than our uncertainties in hadronic and 
nuclear physics. With the assumption ([39]) one finds for the QCD scale [65] 

Aln(A c /M x ) 2vr , 2, 4 , 

l - = ^ d x + 9 d " + 9 ds (40) 

and for the fine structure constant, 

A In a 80a , 43 a , 257 a , 

— = TJa- x dX + 272^* + (41) 



2 We find Aln(m q (Q 2 )/rh q (M%)) = 2/7Alna x ~ (9«x/77r)Aln(A c /M x ) under variation 
of ax, where fh q is the running quark mass and Q 2 — const- A 2 .. 
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Similar values with somewhat different assumptions were found earlier |66|, [T8] . 

For the nucleon mass we include possible strange quark contributions]^ The 
uncertainty in the strangeness content is an indicator of the overall uncertainty 
that may arise due to m s variation. We obtained [I] 

Aln^ = 0.048Aln^ + (0.12±0.1)Aln^, (42) 
A c A c A c 

Aln^ = -0.59Alna + 1.59Aln-^, (43) 
A c Ac' 

where 5n = m n — m p , and thus 

—^t = (0.58 =F 0.08)— + (0.37 T 0M)d s + (-0.65 ± 0.09)d H , (44) 
I ax 

Aln(Gm N ) _ ^ + g ^ Q ^dx_ + , Q ?4 ^ QA1 \ dg + r Q71 ± o.ig)d Hj 
I ax 

(45) 

where the upper or lower signs correspond to the positive or negative signs in 
Eq. (|4*2|) respectively. 

The largest contribution to variations of the proton g-factor g p has been 
argued to arise from the pion loop [28], yielding at first order a dependence on 
the light quark mass of 

Alng p ~ -0.087Alnm/A, 

„ Q M dx_ _ Q7d ^ + QMds (46) 
I ax 

Hence the variations of observables including g p are 

^J^i = (-0.48 ± 0.08)— + (0.59 =F 0.09)d H + (-0.31 ± 0.05)d s 
I ax 

A In v dx 

— _^ = 0.10— - 0.06d H + 0.06d s 

I ax 

AMi? = (1.04 T 0.13)— + (-1.08 ± 0.U)d H + (0.65 T 0.08)c? s 
/ ax 

AMF = (-0.54 ± 0.08)— + (0.65 =F 0.09)d H + (-0.35 ± 0.05)d s . (47) 
/ ax 

We have now expressed the variations accessible to observation in terms of three 
(four) variables: I, dx, dn (and ds), where one parameter may be eliminated 
by normalization. Different unified scenarios will be characterized by different 
relations among these parameters. 

Most data points are upper bounds on a possible variation, with the ex- 
ception of two epochs. First, we consider specifically whether claimed nonzero 



3 In our previous paper [T] we assumed m s /A c = const, here we include the roughly known 
strange contribution to the proton mass. For BBN, the difference in the final dependence 
is less than 3% and hence much lower than the model uncertainty (e.g. for nuclear binding 
energies) . 
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variations of a [28] and \x [37] at redshift 2-3 are compatible with one another, 
since the ratio of their fractional variations is predicted in each scenario. 

Second, we consider whether there is an indication of nonzero variation at 
BBN. For no variation at BBN we obtain x 2 = 17.9 for 3 measured abundances 
( 4 He, D, 7 Li). This discrepancy between theory and observation is exclusively 
due to 7 Li. (Considering only 4 He and D, the value of x 2 is 0.24.) If we wish 
to solve or ameliorate the "lithium problem" by a nonzero variation, we will 
require x 2 l v to be not much larger than unity, taking v = 2 as appropriate 
for one adjustable parameter. If there is no significant range where the three 
abundances have a 2a fit (x 2 ' fv < 4) then we give up the hypothesis that the 7 Li 
problem is solved by coupling variations and instead assume that the observed 
depletion is due to some astrophysical effect. In this case we consider only D 
and 4 He abundances as observational bounds on the size of variations at BBN. 

We will now investigate six different scenarios for the variation of the grand 
unified parameters ax, Mx/Mp, {(p)/Mx and rh/Mx- These will fix the uni- 
fication coefficients dt- For each unified scenario we display the z-dependence 
of the fractional variation (Figs 1-7). Each figure shows the available informa- 
tion from observations of different couplings, interpreted as constraints on the 
variation of a single parameter. These figures are one of the main results of our 
paper. 



Varying a alone Before describing the six different grand unified scenarios, 
we consider a variation of the fine structure constant a alone. Clearly here we 
are unable to account for any nonzero variation in /x or other quantities inde- 
pendent of a. The cosmological history is dominated by the nonzero variation 
of the Ma values at redshifts z ~ 1 to 4. We find that there is almost no 2a 
match of the BBN values {x 2 /v > 3.9): the 2-sigma range is 

3.25% > AlnasBN > 4.06%. (48) 

Hence it seems unlikely that the "lithium problem" can be solved by a variation 
of a alone. If we regard the Li discrepancy as due to systematic or astrophys- 
ical effects we can set a conservative bound on a variation from He and D 
abundances pQ 

- 3.6% > A In dBBN > 1.9%, (49) 

where we imposed that neither the D nor 4 He abundance should deviate by more 
than 2a from observational values. See Fig. [I] for a summary of the bounds in 
this case. 



3.2 Scenario 1: Varying gravitational coupling 

In this scenario we have only du nonvanishing, 

d H = d s = d x = 0, (50) 

therefore 

Ain ^H AinGA - < 5i) 
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Figure 1: Variations for varying a alone. Only observations constraining a 
variation are shown; the BBN fit including 'Li is poor (x 2 /v > 7.8/2) hence we 
also display a conservative bound from 4 He and D abundances neglecting Li. 
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Figure 2: Variations for scenario 2; BBN bounds are 2a bounds. 
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We find that there is no value of AlnGA 2 . for which BBN is consistent with 
the three observed abundances within 2a. The best fit values are x 2 l v > 7.7 
for no variation of mx /Mp at CMB and x 2 l v — 5-9 if the variation of mx /Mp 
has the same size at BBN and CMB. Assuming that the discrepancy in the 7 Li 
abundance is due to some other effect, we find the allowed region of variation of 
G at BBN under which primordial D and 4 He abundance lie within the observed 
range at la (2a), 

- 5% (-13%) < A In GA 2 C < 12% (22%) (52) 
If the variation of mx/Mp has the same size at BBN and CMB one finds 

- 4% (-11%) < A In GA 2 < 10% (16%). (53) 

The bounds on time variation of GA 2 are much weaker than for many other 
varying couplings. This scenario also predicts a vanishing value of rj in Eotvos 
experiments. Thus, to any one of the following scenarios we may add an addi- 
tional nonzero du of similar size to dx, du or ds without changing the results 
significantly. 



3.3 Scenario 2: Varying unified coupling 

In the first GUT scenario without SUSY we consider the case when only dx is 
nonvanishing, 

d x = 1, d H = d s = d M = 0, a x = 1/40. (54) 

Within a supersymmetric theory the same relations will apply except that ax = 
1/24 and the variations of observables are scaled by a factor 24/40 relative to 
A In ax '■ we designate this as Scenario 2S. In both cases we find here 

A In a , . 

— ~ 27. 55 
A In a v ; 

It is then highly unlikely for the nonzero Ma result for variation of a to coexist 
with the determination of [i at redshift around 3 [37], even if the latter is 
interpreted as an upper bound on the absolute size of variation |38j . 

For the BBN fit, we find without SUSY (excluding modifications of the 
baryon fraction r\ due to varying mx) no range of values fitting at la level 
(x 2 / u — 2.3). At 2a the abundances, including 'Li, become consistent for the 
range 

- 5.7 x 10~ 4 < Alnax < -1.7 x 10 -4 (2a). (56) 

If one includes a variation of mx at the time of CMB with the same magnitude 
as at BBN the result remains unchanged (x 2 ' /v > 2.45), with the same 2a 
range. For this scenario we may consider a nonzero variation at BBN, but more 
recent probes must all be viewed as increasingly tight null bounds. 
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3.4 Scenario 3: Varying Fermi scale 

In this scenario we consider the case when the variation arises solely from a 
change in the Higgs expectation value relative to the unified scale, thus only 
dn is nonzero: 

d H = l, d s = d M = d x = 0, a x = 1/40. (57) 
This scenario implies 

^ = "325. (58) 
A In a v ' 

Whether we interpret the determination of fi [37] as a detection or an upper 
bound, any variation in a at large redshift case should be orders of magnitude 
smaller than current observational sensitivity. 

We find for BBN including 7 Li (y = 2) no la range (x 2 /v > 1.95) but 

6 x 1(T 3 < Aln((f))/Mx < 22 x 1(T 3 (2a). (59) 

A variation of mx at the time of CMB with the same magnitude as at BBN 
does not change this result. 

3.5 Scenario 4: Varying Fermi scale and SUSY-breaking scale 

This scenario corresponds to scenario 3, but includes supersymmetry and as- 
sumes that the mass-generating mechanism for SM particles and their super- 
partners gives rise to the same variation: 

d M = d x = 0, d s = d H = l, a x = 1/24. (60) 

We find here 

A In a , . 

— — ^ = -21.5, 61 

A In a v ' 

such that again the claimed nonzero variations in a and \i cannot be compatible 
and the variation in a at redshift 3 must be below current sensitivities. We 
demonstrate this in Fig. HJ where we show for this scenario the bounds on the 
variable dnl = A ln((</>)/Mx) that arise from various observations. Since we 
have only one free variable we can plot all observations simultaneously as a 
function of redshift. Inspection "by eye" permits to judge if a smooth and 
monotonic evolution of dfjl is consistent or not. 

We find for BBN including 7 U(v = 2) no la fit ( X 2 /v > 1.60), while at 2a 

1.25 x 10~ 2 < Aln(4>)/M x < 5.4 x 10~ 2 (2a). (62) 

If one includes a variation of at the time of CMB with the same magnitude 
as at BBN the allowed range becomes slightly restricted (x 2 jv > 1.72), 

1.20 x 10~ 2 < Alm>)/Mx < 4.9 x 10 -2 (2a). (63) 
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3.6 Scenario 5: Varying unified coupling and Fermi scale 

In this scenario we study a combined variation of the unified coupling and the 
Higgs expectation value: 



d M = d s = 0, d x = 1, d H = jd x , a x = 1/40. (64) 

rter 7 cai 
introduced in p] via 



The parameter 7 can be related to the parameter 7 = ^InA^/Mx which was 



2tt 2 

7 = 7 ( 7\ + n7 ) • (65) 



9a x 9 

In pQ we examined the cases 7 = (0, 1, 1.5) which correspond to 7 = (0, 36, 63). 
Here we find that the best BBN fit is reached for 7 ~ 50 with y^/v = 1.45. 
Note that we have the freedom to adjust 7 such that nonzero variations of a 
and [x at redshift ~ 3 are consistent with each other. We have 



A In// 23.2 - 0.657 



(66) 



A In a 0.865 + 0.0027 
We choose for illustration 7 = 42, for which 

A In n = -5.6 A In a (67) 

and the 2a contour for BBN is 

7.5 x 10~ 4 < Alnax < 28 x 10" 4 . (68) 

For a variation of mx at the time of CMB with the same magnitude as at BBN 
the fit becomes worse (x 2 /v > 1.68). However, a 2a fit to BBN is obtained 
over a wide range of < 7 < 26 (negative A In ax ) and 40 < 7 < 00 (positive 
A In ax)- 

Assuming that the apparent 7 Li mismatch at BBN is due to systematic 
astrophysical effects, we may bound ax with only D and 4 He abundances. 
Here we find at la 

- 5.5 x 10~ 4 < Alna x < 1.44 x 10~ 3 (69) 

In Fig. [5] we again plot simultaneously all observations for this scenario. This 
shows that the bound from BBN including 'Li is not consistent with the claimed 
nonzero variations of a and \i for a monotonic evolution over z. 



3.7 Scenario 6: Varying unified coupling and Fermi scale with 
SUSY 

In this scenario we study a combined variation of the unified coupling and the 
Higgs v.e.v. including SUSY, where as in Scenario 4 we tie the variations of the 
superpartner masses and Fermi scale together: 

d M = 0, d x = 1, d s ^d H = jd x , a x = 1/24. (70) 
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Now the relation to 7 is modified as 

7= n^ + ^J (7i) 

One may again adjust 7 to make nonzero variations in a and /i self-consistent. 
With 

Aln^ = 14 - O.287 

A In a 0.52 + 0.0137' 1 ' 

we find that a good fit to BBN is obtained over a large range of 7, ranging from 
7 = 100 to infinity with minimal x 2 l v = 1-45. This shows that the main effect 
in the SUSY model comes from the variation of the Higgs v.e.v. Including a 
variation of at the time of CMB with the same magnitude as at BBN the 
fits gets worse ( X 2 /v > 1.8). A 2a fit can be obtained for < 7 < 28 (for 
negative A In ax at BBN) and for 58 < 7 < 00 (positive Alnax)- 
First, we study the case 7 = 70 for which 

A In n = -3.9A In a (7 = 70) (73) 

and BBN is fit with a 2a range 

5.5 x 10~ 4 < A In a x < 18 x 10" 4 . (74) 

Neglecting 7 Li, we obtain a la bound from BBN 

- 3.5 x 10~ 4 < Alnax < 9.3 x 10" 4 . (75) 
Secondly, we study the case 7 = 25 where 

Aln^ = 8.3Alna (7 = 25), (76) 
and where the 2a contour for BBN is 

- 13 x 10~ 4 < Alnax < -7 x 10" 4 . (77) 

In this second case the Murphy a measurement and BBN point into the same 
direction. The difference between the two values of 7 can be seen from a 
comparison of Figs. and 



4 Epochs and evolution factors 
4.1 Epochs 

In this section we group the information contained in Tables 1-3 and figures 1-7 
into different cosmological epochs. This produces a first quantitative estimate 
of the possible time evolution for the various unified scenarios. The choice 
of epochs is somewhat arbitrary. Two epochs are singled out by events in 
early cosmology, namely the last scattering surface of CMB, and BBN. The 
very recent epoch comprises present day laboratory experiments and the Oklo 
natural reactor, for which a linear interpolation to the present rate of varying 
couplings seems reasonable. We further divide the observations at intermediate 
redshift into three epochs. 
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Figure 7: Variations for scenario 6, 7 = 25; BBN bounds are 2a 



• Epoch 1: Today until Oklo 

Contains Oklo and laboratory measurements. For the laboratory mea- 
surements, we extrapolate the rate of change of the couplings to finite 
changes at the redshift z = 0.14 (t = 1.8 x 10 9 y) of the Oklo event. 

• Epoch 2: 0.2 < z < 0.8 

Contains absorption spectra and isotopic abundance measurements in me- 
teorites. We chose a boundary z = 0.8 since the Murphy dataset [28J has 
relatively few systems around this redshift, making a natural division. 

• Epoch 3: 0.8 < z < 2.4 

Contains several absorption spectra measurements. The end of the Tzana- 
varis dataset [32] sets the cut at z = 2.4. 

• Epoch 4: 2.4 < z < 10 

Contains absorption spectra measurements and bounds on G from neutron 
stars. 

• Epoch 5: CMB, z 1100 

• Epoch 6: BBN, z m 10 10 

4.2 Evolution factors 

We define "evolution factors" l n for epochs n = 1, . . . , 6 by 

AlnG k , n = d k l n . (78) 
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For each unification scenario we will proceed to a quantitative estimate of l n , 
shown in Table The usefulness of considering the evolution factors l n is 
that the unknown (and possibly not monotonic) behaviour of the mechanism 
driving the coupling variations is rolled into a finite number of parameters. 
For a monotonic behaviour they satisfy l n < L whenever z n < z p . The basic 
assumption remains the proportionality AlnGfc(z n ) = dkl{z n ) = dkl n , with 
constant unification coefficients dt independent of the epoch. The normalization 
of l n is arbitrary, and we take for scenarios 2, 5 and 6 

l n = Alna x , n , (79) 
while for scenarios 3 and 4 we take 



l n = Aln((0)/M x ),n- (80) 

For each epoch and scenario, we compute the evolution coefficients l n as a 
weighted average over the measurements in the epoch. The representative red- 
shift z n is the average over the redshifts of observations inside the corresponding 
epoch. It is shown together with the resulting values for l n in table [H This 
table summarizes our results under the assumption of proportionality. 



Rates of time variation in the present epoch 

For Epoch 1 we incorporate the laboratory measurements for rates of varying 
couplings by linear extrapolation in time to the Oklo redshift z\ = 0.14. The 
logarithmic time derivatives may be approximated by linear interpolation 

— = d t lnG k ~ — -, (81) 

Lrk to — El 

where t\ = 1.8 x 10 9 y is the time corresponding to the redshift z\ = 0.14. 



Method of averaging 

We evaluate the weighted average using all values listed in table [TJ This pro- 
cedure may be quite problematic, since sometimes different observations are in 
manifest contradiction. We take the attitude that, given the possible presence 
of systematic effects both in spectroscopic determinations of nonzero coupling 
variations and in the primordial 7 Li abundance, a viable model need not fit all 
data points. However, even if any given nonzero claimed variation is actually 
due to systematic error, we still expect the size of the error to be comparable 
to the size of the claimed variation. Thus, such claims are most conservatively 
interpreted as bounds on the absolute magnitude of variation. The surviving 
nonzero variation(s), in addition to the null bounds at other epochs, define a set 
of evolution factors which must be satisfied by any explicit model of evolution. 

For some scenarios we therefore also evaluate the evolution factors that 
are obtained by considering that some of the claimed observations of nonzero 
variation may instead be due to an underestimated systematic error. These 
alternative evolution factors are given in square brackets, corresponding to the 
following replacements: 
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Epoch 


1 


2 


3 


4 


5 


6 




0.14 


0.53 


1.6 


3.8 


10 3 


10 10 


Scenario 


h x 10* 


l 2 x 10 e 


h x 10 b 


h x 10 s 


h x 10 4 


k x 10 3 


a only 


-0.01 ± 0.06 


-1.1 ± 1.0 


-0.26 ±0.10 


-0.85 ±0.37 


-150 ±350 


5 ±34 


2 


-0.1 ±0.1 


0.04 ±0.03 


-0.15 ±0.08 


0.10 ±0.03 


0.9 ± 14 


-0.37 ±0.20 


3 


4.1 ±4.8 


-1.5 ± 1.2 


0.42 ± 3.3 


-3.6 ±0.9 


69 ± 920 


14 ±8 


4 


3.9 ±8.5 


-3.4 ± 2.7 


-8.4 ±5.1 


-8.7 ±2.1 


31 ± 450 


33 ±21 


5, 


-0.02 ±0.18 


-0.24 ±0.18 


-0.25 ±0.10 


-0.61 ±0.13 


0.6 ±8.6 


1.7 ± 1.1 


(7 = 42) 












[0.4 ± 1.0] 


6, 


-0.02 ±0.12 


-0.10 ±0.07 


-0.17 ±0.07 


-0.44 ±0.10 


0.3 ± 5.0 


1.2 ±0.6 


(7 = 70) 












[0.3 ± 0.6] 


6, 


-0.12 ±0.18 


0.04 ±0.12 


-0.30 ±0.11 


0.29 ± 0.08 


0.7 ±10 


-1±0.3 


(7 = 25) 








[-0.43 ±0.28] 







Table 4: Redshifts and evolution factors for each epoch, for the scenarios defined 
in section [3j In the first row the values of l n give the fractional variation of a; 
in Scenarios 2, 5 and 6 that of ax] and in 3 and 4 that of ((f)) /Mx- Values in 
brackets give, for BBN (Iq) the evolution factors neglecting 7 Li; or for I4, the 
evolution factor with the A/i/fi value of [37] substituted by that of [55] , 



Scenario 5, 7 = 42: Neglecting Li-abundance at BBN 
Scenario 6, 7 = 70: Neglecting 7 Li-abundance at BBN 

Scenario 6, 7 = 25: Replacing the fi measurements of [37] by the conservative 
upper bound of [38] . 

In the case where a alone varies, since the fit including 7 Li is poor we calculate a 
2a range using observational central values and errors of D and 4 He abundances 
given in pQ. 

4.3 Monotonic evolution with unification 

Here we briefly summarize whether the unified scenarios we consider can be con- 
sistent with a monotonic evolution of the single underlying varying parameter, 
based on the evolution factors l{ found in Tabled! 

Varying a only 

Although variation of a alone does not help to account for deviation of BBN 
abundances from standard theory, or for any nonzero variation of fi, the cosmic 
history is interesting due to the significant nonzero value in Epochs 3 and 4. The 
Oklo bound in Epoch 1 restricts the present time variation to 3.7 x 10~ 17 y _1 
(assuming no acceleration of dta). 

Scenario 2 

Scenario 2 favours a negative variation of ax at BBN, and a negative variation 
may also fit the Ma results. However, the Reinhold fj> measurement indicates a 
positive, but much smaller, variation. We keep the R^ results, which dominate 
the weighted average due to their small error on A max, to obtain Z4. The 
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ratio A In /i/A In a = 27 makes this scenario unlikely to fit the reported signal 
of nonzero Aa. 

Scenario 3 

In scenario 3 a positive variation of {4>)/Mx is favoured by BBN. The high ratio 
Aln/i/Alna ~ —325 makes the bounds obtained on a variation of \x strongly 
inconsistent with the claimed size of variation of a. We keep the Reinhold et 
al. values to obtain I4, which again dominate the results. 

Scenario 4 

In this scenario, the ratio A In ///A In a = —22 is again large and makes any 
observation of significant nonzero A In a unlikely. Both the Ma and the R/i 
measurements point in opposite direction to BBN; however the two spectro- 
scopic observations are also inconsistent with each other, within the scenario. 
Again, we keep the R/i results which dominate the determination of I4 due to 
the small error. 

Scenario 5, 7 = 42 

In this scenario the variation of ax favoured by BBN is positive (Iq = (1.7 ± 
1) x 10 3 ), however both nonzero variations from spectroscopic data Ma and 
R/U require negative variations. With Aln/i/Alna = —6 the spectroscopic 
measurements appear consistent with each other. Hence one would require 
some non-monotonic evolution to fit nonzero variations both at BBN and at 
moderate z. Thus in TableHJwe have also evaluated Iq using only the constraints 
given by D and 4 He (in brackets). 

Scenario 6, 7 = 70 

As in the preceding scenario, BBN favours a positive variation in ax, but Ma 
and R/i favour negative. Again, Fig. [6] may suggest a non-monotonic evolution. 
Fitting to BBN including 7 Li we would obtain Z 6 = (1.2 ± 0.6) x 10~ 3 ; table H 
also displays in brackets the value of Iq obtained from D and 4 He bounds only. 

Scenario 6, 7 = 25 

In this scenario, both BBN and the Ma signal favour a negative variation of ax, 
whereas the R/i observations point towards a positive variation. Following the 
argument of Wendt et al. [38J , we substitute the R/i value by the null constraint 
|A/i//x| < 2.5 x 10~ 5 [38j to obtain the bracketed value of Z4 in Table HI In this 
scenario the evolution factors show a crossover from negligible variation at low 
redshift, to strong and monotonically increasing negative variation at z ~ 2. 

4.4 Tension between the 7 Li problem and variation of // 

Measurements of the primordial 7 Li abundance show that the BBN abundance 
needs to decrease below the standard value to fit the observations, whereas the 
Reinhold /i measurement indicates /i to increase at z ~ 3. We find that for all 
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our unification scenarios the sign of the dependence on the fundamental param- 
eter is the same for \i and Li. Moreover, the coefficients of this dependence 
are nearly identical up to a common factor; hence the induced variations for \x 
and 7 Li point in the same direction, in contradiction to the tendency inferred 
from the observations. For example, for scenario 5 we find 

A In \t = (23.2 - O.657) A In a x , 
Aln 7 Li= (1692 -497) A In a x - (82) 

These expressions change sign at 7 = 35.7 and 34.5, respectively. For a mono- 
tonic evolution, there is no possibility to have both a significant variation of \i 
and a variation of opposite sign in the Li abundance. (In the regime 7 ~ 35 
there is no 2a fit to BBN.) A similar result can be found for scenario 6 (including 
the SUSY partner mass dependence, which shows the same sort of degeneracy). 
Note that scenario 2 and 3 are just limiting cases of scenarios 5 and 6. 

The main reason for this behaviour is that variations of 7 Li and \x are domi- 
nated by the variations of m/A c and m e /A c , respectively, with the same sign of 
prefactor. This degeneracy can be broken if m e varies differently from the quark 
masses, a possibility that we do not consider in this paper. For our scenarios 
with constant m/m e , the conflict between a monotonic time evolution and the 
\i- and Li-observations is reflected in the opposite signs of Z4 and Iq. 

This observational tension for monotonic behaviour is clearly depicted in 
Fig. El where we plot simultaneously the averaged observational values of evo- 
lution factors /j/ ln(l + Zj), normalized to Z4/ ln(l + 2:4). For Scenario 6, 7 = 25, 
we also display the result obtained by substituting the Wendt et al. value of \x 
variation for that of [37J. The factor ln(l + Zi) is introduced as a convenient 
normalization to avoid compressing the scale of variations excessively in recent 
epochs@ For the purpose of a quick inspection we have omitted the error bars, 
which are of course necessary for a quantitative interpretation. 

4.5 Special values of 7 

In Scenarios 5 and 6 there is a value of 7 for which Aln(</>)/A c vanishes. For 
these values, Standard Model physics undergoes an overall multiplicative shift 
of energy scale under variation of ax, up to variations of perturbative, dimen- 
sionless couplings: specifically the Yukawa couplings (whose variation we have 
generally neglected) and a. The significant observable effects arising from vari- 
ation of SM fermion masses relative to A c , which dominate in most unified 
scenarios, are largely absent, and the low-energy phenomenology is very similar 
to the case of varying a only. In particular the 7 Li problem at BBN is not 
addressed and the variation of \x is smaller than that of a. 

The required values are 7 = 2ir/7ax — 36 in the case without SUSY (ax — 
1/40); or 7 = 27r/3ax ~ 50 with SUSY (ax — 1/24) when the superpartner 
masses vary with the Fermi scale, ds = du = jdx- From a low-energy point 

4 In quintessence-like theories, if the scalar field contributes a constant fraction of the total 
energy density of the Universe, as in so-called "tracker" models, the evolution of the field is 
typically also proportional to ln(l + z). This is an additional motivation for our normalization. 
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Figure 8: Normalised evolution factors U/I4, for each scenario, where li = 
k/lnil + Zi). 

of view these values appear as fine-tuning, however it is conceivable that they 
would arise from some specific mechanisms of electroweak symmetry-breaking 
or SUSY-breaking. 

5 Summary / Conclusions 

Within Grand Unified Theories, different measurements of the variation of fun- 
damental constants can be consistently reduced to a variation of a few "unifi- 
cation parameters", namely the unification scale Mx/Mp, gauge coupling ax, 
the Fermi scale {4>)/Mx and SUSY-breaking masses rh/Mx- We define various 
GUT-scenarios for varying couplings by the assumption of proportionality of 
fractional variations of the unification parameters. 

Assuming that couplings really vary, this is a way of excluding such GUT 
scenarios by demanding consistency of the implied variations. The assumption 
of proportionality permits us to project all observations into constraints on a 
common evolution factor l{z) for each scenario. We show that different GUT 
scenarios yield different time evolutions of l{z) assuming that certain claimed 
measurements of varying constants are correct. We confirm that "simple" mod- 
els which have only one fundamental parameter varying {ax or {(f)) /Mx) result 
in inconsistent variations. However, combined variations of these two parame- 
ters, as described in scenarios 5 and 6, lead to results more consistent with the 
possible quintessence-induced time variations of fundamental couplings which 
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we investigate in [7J. 

Specifically, one may ask whether the claimed observations of variations in 
a. [2S] and (i |37j are mutually consistent, and whether they are consistent with 
an explanation of the apparent primordial 7 Li-depletion by varying couplings. 
Within a hypothesis of constant Yukawa couplings, which results in identical 
fractional variations of all quark and lepton masses, we investigated arbitrary 
variations of ax, ((f)) /Mx and Mx/Mp. For scenarios with supersymmetry 
we also assumed that the SUSY-breaking masses vary proportional to ((f)), but 
the effect of such a variation only appears at higher order and is probably not 
crucial. 

We have not found a scenario with a monotonic time evolution l(z) that 
makes all three signals or hints of variation mutually consistent. A monotonic 
evolution requires either to discount one of the "signals" by substantially in- 
creasing its uncertainty, or to alter our assumptions by including additional 
time variation of some Yukawa couplings. 

Our investigation shows how the variations of different couplings in the 
Standard Model may be compared. If the observational situation becomes 
clearer and at least one nonzero time variation is established, such methods 
may be used for new tests of the idea of grand unification. 

Note added 

Shortly before the completion of this paper a new determination of the variation 
of fj, appeared [67] reporting a reanalysis of spectra from the same two H2 
absorption systems as [37J, and adding one additional system at z ~ 2.8. The 
results of the new analysis are not consistent with the previous claim indicating 
a nonzero variation, either considering all three systems or the two previously 
considered. The stringent null bound of the new analysis, Afi/fj, = (2.6 ±3.0) x 
10 -6 , would disfavour all scenarios except those where the fractional variation 
of fi was of the same order as or smaller than that of a. This would require 
us to approach the "special", apparently fine-tuned values of 7 discussed in 
Section 14.5^ for which /x variation (and any deviation from the standard 7 Li 
abundance at BBN) are suppressed. 
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Appendix A Effect of "varying constants" at CMB 
and 77 

In our previous work on BBN we used the WMAP determination of the baryon 
number density parameter r\ = ns/n 7 directly to reduce by one the number of 
unknown parameters. However, we should also consider the effect of possible 
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variations of Gk at the epoch of CMB decoupling. This question has distinct 
aspects: first, can the CMB alone or combined with various other cosmological 
observations give useful bounds on the values of fundamental parameters at this 
epoch? Second, how do the possible variations affect the determination of 77? 

It would not be appropriate to give an extended discussion of CMB bounds 
on fundamental variations here; the subject has already been treated [25] at 
length. Bounds tend to depend strongly on the values taken by cosmolog- 
ical parameters which are not at present well known through independent 
measurements: in other words there is considerable degeneracy. Fundamen- 
tal parameters affecting the CMB are the proton and electron masses, the 
gravitational constant and the fine structure constant, as well as the mass of 
any dark matter particle present. In Planck units, these reduce to the par- 
ticle masses and a. The relevant cosmological parameters are the amplitude, 
spectral index (and possible running, etc.) of primordial perturbations; the 
baryon, dark matter and dark energy (cosmological constant, etc.) densities 
normalised to the critical density; the Hubble constant; and the reionization 
optical depth. Of these, the baryon density ft^h 2 will vary linearly with the 
proton mass in Planck units, for a fixed baryon-to-photon ratio 77. Conversely, 
given a measurement of f^/i 2 , the correct value of 77 varies inversely with the 
proton mass. The conversion factor between Q^h 2 and 7710 = 10 10 ?? is then 
273.9{m p VG) l0 {m p VG)- 1 ~ 273.9(1 - A ln(m N /M P ) l 

cmb)> where we approxi- 
mate the proton and neutron masses by their average mjv- 

If, therefore, we allow the proton mass (or the gravitational constant, in 
QCD units) to vary arbitrarily at the CMB epoch, 77 is undetermined by WMAP 
and we must consider it as an extra free parameter or try to impose independent 
cosmological bounds. However, we impose that the size of variations away 
from the present value of m p /Mp is a monotonically decreasing function of 
time: thus A ln(m7v/Mp)| CMB < A ln(mAr/Mp)| BBN . Hence we would have 
a self-consistent treatment of this parameter if the secondary discrepancies in 
primordial abundances due to an incorrectly estimated 77 were smaller than the 
primary effect of varying /Mp at BBN. The relevant results of our previous 
analysis 



dln(D/H,y p , 7 Li/H) 
<91n(mAr/M P )| BBN 
ain(D/H, Y p , 7 Li/H) dln(D/H, Y p , 7 Li/H) 



<91n77 d\u{rriK[/Mp) 



CMB 



(1.88,0.72,-1.14), 
(-1.6,0.04,2.1) (83) 



are derived in QCD units where the strong coupling scale A c is constant, and 
where we neglect small contributions to the nucleon mass tun and take it pro- 
portional to A c . The first relation, derived at a fixed value of 77 = 6.1 x 10 -10 
(WMAP3 [68] jH led to the bound -0.095 < A ln(mAr/M P )| BBN < +0.05, 
where the main sensitivity to this variation is due to helium-4 (Y p ). Since 
this abundance is insensitive to changes in 7/, we postulate also that —0.095 < 
Aln(m7v/Mp)| CM B < 0.05. 



1 Updating to WMAP5 values does not lead to any significant change 
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The resulting errors in the (standard) BBN abundances due to a possibly 
misestimated r\ are then 

5 ln(D/H, Y p , 7 Li/H) = ({-0.15, 0.08}, {-0.002, 0.004}, {-0.10, +0.20}) (84) 

to be compared with observational errors of 

<td/(D/H) ~ 0.4/2.6 ~ 0.15 
vme/Y p ~ 0.009/0.25 ~ 0.04 

CT 7Li /( 7 Li/H) ~ 0.5/4.5 ~ 0.1, (85) 

where we take the standard BBN 7 Li abundance 4.5 x 10~ 10 as central value. 
Hence the variation of mpf/M-p at the CMB epoch and consequent rescaling 
of rj may in principle have significant consequences for deuterium and lithium 
abundances in BBN. It may be appropriate to take A ln(mjv/-Wp)|CMB as an in- 
dependent variable in the analysis of BBN variations. The maximum effect due 
to rescaling of r\ would occur when Aln(mAr/Mp)| CMB = A ln(mAr/Mp)| BBN , 
giving a total sensitivity of 

^gMS _ (3.48, 0.68, -3.24). (86) 
dln(mAr/Mp)| BBNiCMB 

Appendix B The 8 Be resonance 

The 7 Be+n — > 7 Li+p reaction is the main channel for destruction of 7 Be during 
BBN. If this reaction was not present, the final 7 Li abundance predicted by 
standard BBN would be considerably higher: 

7 Li/H = 4.5 x 10~ 10 ~ 14 x 10~ 10 . 

The high cross section of this reaction is due to a strong 8 Be resonance which 
sits at about the energy of both 7 Be+n and 7 Li+p [69J. For the reaction to 
continue to operate efficiently, it is important that the resonance remains near 
these 7 Be / 7 Li energy levels. We will argue here that, given the size of coupling 
variations relevant for our paper, this is indeed the case. 

In [1] we estimated the dependence of nuclear binding energies on the pion 
mass by 

|^ = MAi - l)^°r * -0.13/i(A - 1), (87) 

taking r ~ —8. The constants /» are expected to be of order unity, but will 
differ between light nuclei due to peculiarities of the shell structure. Our nor- 
malization corresponds to /d = 1- We are then concerned with the relative 
changes of the Be and Li binding energies and the energy of the 8 Be reso- 
nance, whose dependence we will estimate in an analogous way with a constant 
of proportionality /§. Then 

AB 7Be = -9.1 MeV x 6/ 7Bo Alum, 
A£ 8Be . = -9.1 MeV x 7/g Alnm, 
AB 7Li = -9.1 MeV x 6/ 7Li Alnm, 
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recalling that to^ oc to 1 / 2 . In [70] the sum of the neutron and proton widths of 
the 8 Be resonance is given as approximately 1.6 MeV thus for the 8 Be-destroying 
reaction to remain effective we require at least 

|9.1MeV(7/g -6/7) A In m| < 1.6 MeV, (89) 

where fj may correspond to either 7 Be or Li. If we take all fi = 1 this condition 
becomes | A In rri | < 0.18, easily satisfied by the range of variations that we 
consider (A In to was bounded at about 1.5%). However, this would imply a 
substantial cancellation between the variations of A = 7 and ^4 = 8 states, 
which may not occur for the true values of fi. There may be less cancellation, 
for example if /g = 2, = 0.5 we obtain |Alnm| < 0.016, which is still fulfilled 
in the unified scenarios we consider where the variation of to is around 1%. 
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